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■ 1 Introduction 



The present status of the relativistic direct interaction theory fl], ||, |3|, £| enables us 
to consider it as a natural basis for a consistent description of relativistic effects in 
various physical systems on different levels: classical, quantum, or statistical ones. 
The relativistic statistical mechanics of interacting particle systems is nowadays at 
an early stage of its development, although the classical partition function of a rela- 
tivistic ideal gas was calculated by Jiittner in 1911. One cause of this was indicated 
in report [Q. Ter Haar and Wergeland wrote: "At extremely high temperatures 
relativistic effects may, of course, be important. Then, however, matter behaves as 
a mixture of ideal gases and this limiting case poses no problem. By and large, 
a relativistic theory of heat seems, therefore, to be of little practical importance" . 
But the 'practical importance', being a very nonsmooth function on historical time, 
cannot be considered as the main reason for theoretical investigation. 
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Among various approaches to the relativistic direct interaction theory the single- 
time Lagrangian formalism [[!], ^ |7J seems to be the most convenient in the consider- 
ation of the general problem of relativistic dynamics, as well as in the investigation 
of various approximations. This formalism has been extended to an arbitrary form 
of relativistic dynamics defined geometrically by means of space-like foliations of the 
Minkowski space ||, ||. Transition from the classical Lagrangian to the Hamilto- 
nian description allows one to consider the relativistic effects in the statistical and 
quantum mechanical properties of the particle systems. 

The present paper is concerned with relativistic models of the equilibrium sta- 
tistical mechanics. In Sec. 2, a brief introduction into the concept of the form of 
relativistic dynamics is presented. Investigation of a classical and quantum relativis- 
tic ideal gas by means of the front form of dynamics is outlined in Sec. 3. Relativistic 
generalization of a one-dimensional gas of hard spheres is obtained in Sec. 4. This 
gives an example of an exactly solvable model in the relativistic statistical mechan- 
ics with a nontrivial particle interaction. Sec. 5 is devoted to the investigation of 
the first quasirelativistic (post-Newtonian) approximation. In the instant form of 
dynamics the weakly-relativistic corrections to the thermodynamical functions of 
the dilute gas with short-range interactions are studied on the basis of the general 
structure, to the order c -2 , of the approximately relativistic Hamiltonian. As an 
example, the relativistic correction to the Van der Waals equation to the order c~ 2 
is obtained. 



2 Forms of relativistic dynamics in the Lagran- 
gian description of an interacting particle sys- 
tem 

Let us consider a dynamical system consisting of N interacting point particles. It 
is convenient to describe the evolution of this system in the (n + l)-dimensional 
Minkowski space M n+1 with coordinates x^, \i = 0,1,..., n. We use the metric 
WV/wW = diag(l, —1, • • • , — 1 ). In applications we put n = 3 or n = 1. The motion of 

n 

the particles is described by the world lines 

7a :R->M n+ i, r ^<( Ta ); a = l,...,N, (2.1) 

being time-like one-dimensional unbounded submanifolds of the Minkowski space. 

The relativistic freedom in the simultaneity definition makes possible different 
three-dimensional descriptions of the relativistic particle motions. According to 
Dirac [lCj, p| , they are called the forms of relativistic dynamics. Within the frame- 



work of the single-time Lagrangian or Hamiltonian mechanics this concept may be 
introduced in the following way || [7| . Let us consider foliation E of the Minkowski 
space M n+1 by the hypersurfaces 

t = a(x), t el, (2.2) 
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with the next property: every hypersurface S 4 = {x <G M n+1 | a(x) = t} must 
intersect the world lines 7 a of all the particles in one and only one point 

x a {t)= la ^ t . (2.3) 

This allows us to consider t as an evolution parameter of the system [Hp. In the 
Poincare-invariant theory, when we consider only time-like world lines, the hyper- 
surfaces ( |2.2j ) must be space-like or isotropic: 

V^aWa) > 0, (2.4) 



where <9 M = d/dx^. Then, we have d°a > 0, and the hypersurface equation ( |2.2j ) has 
the solution x° = (p(t, x), where x = (x l ), i = l,...,n. Therefore, the constraint 
x a (t) G S t enables us to determine the zeroth component of x a (t) in terms of t and 
x l a (t). The parametric equations ( |2.1| ) of the world lines of the particles in the given 
form of dynamics have the form: 

a; = <p(t,x a (t))=<p a , x l = x\(t). (2.5) 

The evolution of the system is determined by nN functions t \— > x l a (t). They 
may be considered as representatives for the sections s : R — > F, i m (t,x l a (t)) of 
the trivial fibre bundle 7r : F — > R with niV-dimensional fibre space E = M. nN . The 
latter constitutes the configuration space of our system. 

Three Dirac's forms of relativistic dynamics correspond to the following hyper- 
surfaces (|2.2|) : x° = ct (instant form), x° + x n = ct (front form), and r] liu x^x v = c 2 t 2 
(point form). Other examples for the case n = 3 may be found in 0. 

In the relativistic Lagrangian mechanics the Lagrangian function L : J°°7i — > 
R is defined on the infinite order jet space of the fibre bundle n : F — > R with 
the standard coordinates x 1 ^ ||12|| . The values of these coordinates for the section 
s : t t— > (t,x l a (t)) belonging to the corresponding equivalence class from J°°n are 
xf\t) = d s x{{t)/dt s = D s xi, s = 0, 1, 2, . . .. 

The free-particle system is determined in any form of dynamics by the Lagrangian 
Lf : J l n — > R depending on the first derivatives: 



L f = J2 macy/ip^t, x)) 2 - v* a = xf\ (2.6) 

a 

m a being a rest mass of the particle. In the front form of dynamics it reads: 



L f = 5> aC yi - (v 2 al + ■■■ + ^ (n _ 1} )/c 2 - 2v an /c. (2.7) 

a 

In the general case the Poincare-invariance conditions forbid the existence of 
interaction Lagrangians which are defined on the jet-space J r ir with some finite r 
(for example, with r = 1). This leads to serious difficulties in physical interpretation 
of the formalism, and, in fact, makes it impossible to obtain a closed form of the 
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corresponding Hamiltonian functions which is necessary for the development of the 
classical statistical mechanics []I|, [7|, [8J. This fact is the Lagrangian counterpart of 
the famous no-interaction theorem in the Hamiltonian relativistic mechanics |13| . 



There are at least two possibilities to avoid this difficulty. The first is offered 
by the front form of dynamics in the two-dimensional Minkowski space. In this 
case there exists a wide class of interaction Lagrangians depending on the first 



order derivatives |T^] . This allows us to obtain corresponding Hamiltonian functions 
by the standard Legendre transformation that preserves the physical meaning of 
the position canonical coordinates. The second consists in the consideration of 
the approximation in c~ 2 Jl|, ||, [15[]. At least in the first order approximation we 



may obtain the usual Lagrangian functions depending on the first order derivatives 
IJ. Examples of using these possibilities in the investigation of relativistic effects 
in models of the equilibrium statistical mechanics will be considered in the next 
sections. 



3 Relativistic ideal gas within the framework of 
the front form of dynamics 

As an illustrative example we shall consider here a description of the relativistic 
ideal gas by means of the front form of dynamics. The Hamiltonian function of the 
free iV-particle system has the form: 

N 

H = J2Ha(*a,Pa) (3.1) 

a=l 

with 

g a = ffo(p) = c m Y + Pl . (3.2) 

This Hamiltonian can be obtained from the Lagrangian ( |2.7| ) by the usual Legendre 
application. 

The n-dimensional coordinates x a cover the given region Q C W 1 with vol(fi) = 
V, and p a belongs to the region 77 = {p a = (p al , . . .,p an ) G R n \ p an > 0}.As it is 
well known, all the information about thermodynamical properties of the system is 
contained in the expression for the canonical partition function fT6| : 



h nN N\ 



J e-? H l[d n x a d n p a , (3.3) 



where, as usual, (3 = (kT) 1 and the integration is performed over the phase space 
P. In our case P = Q N x II N . Inserting ( |3.1| ) into Eq. ( |3.3| ) we obtain 



7jn ~ h^W\ zN ' (3 " 4) 
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where 



z = J e~ l3Ha d n x a d n Pa . (3.5) 
Substitution of the front form Hamiltonian (|3.2| ) immediately yields 



v [ dPn e-^ + ^) n / d Pk 

n k—1 „ 



g 2p n ^fc _ 



(3.6) 



Performing integration over p^. and putting p n = mca we have 

z = Vmc (n " V2 Jdaa^^l 

o 

Using the integral representation of the Hankel function K n (x) = He^ffn'i 



(3.7) 



IX 



K n {x) = ^J daa n - x e-^ a +^, 
o 



x > 0, 



we obtain 



z = 2Vmc 



( 27rm 



(n-l)/2 



K {n+1)/2 {f3mc 2 ). 



When n = 3, Eqs. ( |3.4j ) and ( |3.9| ) give well-known Jiittner's result: 



ZiV " iV! 



47rm 2 cV 
(3h 3 



-i N 



K 2 (pmc 2 



which is usually derived by means of the instant form Hamiltonian: 

H a = c\/m 2 c 2 + p^. 
Using the asymptotical expansion 



KJx) 




7T / 4n 2 - 1 

e _!t H h . 

2x V 8x 



x — > oo, 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



(3.12) 



we can obtain a weakly- relativistic correction to the nonrelativistic result: 



7 qo) p Nm ? 1 1 , < n + 2 ) 
Zat - Z N e 1 + — — — 

5 8pmc z 



N 



where 



7 (0)id 



VA 



-nN 



A = 



(3h 2 
2nm 



(3.13) 



(3.14) 



The necessity for (rather trivial) renormalization of the nonrelativistic partition 
function zffl — > zffl e~^ Nmc2 follows from the presence of the rest energy mc 2 in the 



relativistic Hamiltonians ( |3.2| ) or (|3.11| ). 

The obtained expressions can be useful for the treatment of the quantum rela- 
tivistic ideal gas with generalized statistics of arbitrary order q. In this case it is 
convenient to consider the grand partition function 



9 = Tre-^ H ^ N) = e~ pn , 



(3.15) 



where H and N represent the operators of energy and of particle number, respectively, 
and n denotes chemical potential fl6|] . 

Assuming that the maximum value of the occupation number in a given energy 
state can be q, any integer greater than 1, and acting in a manner quite similar to 
the nonrelativistic case (cf., for example, |T7|]), we obtain: 



gv_ 

' f3h n 



d n p\n 



-0(ff+l)(flb(p)-/O 



(3.16) 



Here g denotes degeneracy of the energy state with a given value of momentum p. 
For structureless particles with the spin a we have g = 2a + 1. 

Of course, care is needed in the replacement of the summation by the integration 
over the momentum states, when q ~ N, and, especially, if q — > oo (that corresponds 
to the Bose-Einstein statistics). But this problem, connected with the Bose conden- 
sation, is essential at a low temperature, << mc 2 , when relativistic correction 
seems to be negligible. Then the standard thermodynamical relations lead to the 
following expression for the pressure: 

P = gfi- x hT n J d n p [In (1 - A' ?+1 e- /3{9+1)ffo(p) ) - In (l - Ae _/,Ho(p) )] , (3.17) 



where 



is the fugacity. The number density 



A = e^ 



is given by 
1 



gh~ n / d n p 



1 _ N _ dP 
v V dfi 



^ e -/3#o(p) _x g+1 e' /3(9+1) ' ffo(p) 
-(g + 1) 



(3.18) 



(3.19) 



1 _ \ e -0H o (p) ^ I 1 _ \q+l e -P{q+l)H {v) 

and the inner energy density u = U/V is found to be 



(3.20) 



u = g hr n \ d n P H (p 



\ e -PH ( P ) _x 9+1 e" /3(9+1) - ffo(p) 

-(q + i)- ' 



1 - Ae-^o(p) 



1 _ \q+l e -(3(q+l)H (p) 



(3.21) 
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Expanding the exponential functions into the series and using the denotation 

<j>{(3) = gh~ n z/V, (3.22) 
where z is defined by Eq. (|3.5|), we have 



f3P = k~ l [\ k (f){f3k) - A fc(9+1) 0(/3(g + 1)A;)] (3.23) 

k=l 

and 

1 oo 

- = J2 [avw -( q + i)\ k ^mq + m] ■ (3.24) 

k=l 

Keeping only the terms up to the second order in the density 1/v, we obtain an 
equation of state of the form: 

where 

"« = {t''>i. (X26) 

Similarly, for the energy density with the same accuracy we have 



Making use of the explicit expression for the function from ( |3.9|) and ( |3.22| ), we 
may find a relativistic equation of state at a high temperature up to the first order 
in the degeneracy parameter 

A n 

5=—. (3.28) 

gv 



The equation is found to be 

v 

where 



(3.29) 



Because F(x) —>■ 1 as x — > oo, Eq. ( |3.29| ) agrees with the known nonrelativistic 
results [fjjjfl . In a similar manner we may write relativistic virial expansions for other 
thermodynamical functions of the quantum ideal gas, especially for the inner energy 
and specific heat. 
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4 Relativistic one-dimensional model of the hard 
spheres gas 

In the two-dimensional space-time M 2 , the front form of dynamics serves a wide 
class of nontrivial interaction Hamiltonians in terms of the covariant canonical co- 
ordinates. 

The general form of a Hamiltonian function for a system of N identical particles 
on the line I = {xGlR|0<x< A} described by the canonical coordinates x a and 
momenta p a within the framework of the front form of dynamics is given by |T0] : 

H = J^HoiPa) + ^^{Pa+Pb)V{r ab p a ,r ab p b ). (4.1) 

a a<b 

Here 

H (p) = Up+—), (4.2) 



2 V P 

r ab = x a — x b , and V is an arbitrary function on the indicated arguments. In this 

section we put c = 1. From the definition of the front form of dynamics it follows 
that momenta p a belong to the positive semiaxis M + = {x 6 R | x > 0}: 

Pa > 0. (4.3) 

For the convenience of the comparison with the corresponding nonrelativistic calcu- 
lations we assume that function V has the form: 

V = V(|?«6|), (4.4) 

where 

q a b = r ab {p a + Pb)v ( — + — ), (4.5) 

\Pb Pa) 

and v : R+ — > M + is some function which will be defined later. 
Let us again consider the canonical partition function 



j^yj J e-? H l[dx a dp a , (4.6) 



over the phase space P = I x R^. 

Since we are interested in the relativistic generalization of the hard spheres 
model, we can choose interaction function (14.40 in the form 



V(x) = | 


o, 


x > a 




x < a. 


oo, 


PV(\q ab \) = 




\q a b\ > o 






\q a b\ < o"- 



(4.7) 

We, therefore, find 

(4i 



S 



If we choose function v(x) in such a way that for any a, b, c 

if \q ab \ > a and |g 6c | > er, then |g ac | > a, (4.9) 
the partition function ( |4.6| ) can be written in the form: 

a a 

where 

d/x(p) = e- m){p) dp (4.11) 
and coordinates x a must satisfy the restrictions 

M>^— uf^ + ^V tt^^-^r. (4-12) 

Pa+Pb \Pb Pa J V{X) 

Conditions ( |4.9|) lead to the inequality 

X - < 44 < — ■ (4-13) 

4 - u(2) ~ 4 V ; 

In the following we shall choose the simplest solution of (|4.13|) which has the form 

u(x) = C(x + 2), (4.14) 

such that 

g a6 = C' 1 ^ (- + - ) . (4.15) 



Because in the nonrelativistic limit p a — > m, the demand that q a b — ► r a b in this limit 
fixes the value of the constant C: 

C=j. (4.16) 
Therefore, partition function (|4.10Q may be rewritten in the form: 

z " = / (f[ MrMs) n e - *. - ^ - . (4.17) 

The appearance of Heaviside 6 functions is a consequence of conditions ( [4.12 ). Let 
us perform in ( |4.17| ) a change of the variables (x a ,p a ) — > (y a ,p a ), such that the 
arguments of 6 functions become y a +i — y a - This gives 

ma 

Va = %a —fa, (4.18) 

where 

pi = o, & = — + —, (4.19) 

Pi P2 
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and 



i £i i i 

<p a = — + 2 \ - + — 

Pi ~ Pb Pa 



a>2. 



(4.20) 



As a consequence, we obtain 



Z N (A,(3) 



1 r N 

a=l 



ma \ w 



(4.21) 



Before performing an integration over momenta we consider the Laplace transfor- 
mation of ( |4.21|) 



Zn(s,(3) = J dAe~ sA Z N (A,P) 
o 

i r N 

= h N s N+l / n d KPa)ew(-S^- 



(4.22) 



Using expressions ( 4.11|) and ( 4.19 ), ([4.20 ) we find 

1 



Z N (s,/3) 



h N s N+l 



z(/3,M 2 ) 2 z((3,M 1 ) 



N-2 



where 



and 



Ml = m 2 + 2s- 



ma 



Ml = m 2 + s 



ma 

T'' 



DC 

z(f3,m) = J dpexp 



m 



p 



2mK 1 ((3m) 



Next, we consider the grand partition function 



(4.23) 
(4.24) 

(4.25) 
(4.26) 



A' 



The summation over N is performed immediately giving 

z 2 ((3M 2 ) h 



Z(J3,8,X) 



(4.27) 



z 2 ((3M x ) sh - Xz((3M 1 ) ' 
The asymptotic behaviour of function Z(j3, A) is determined by the singularity points 



of function ( 4.27|) which lie on the real axis for variable s. In our case there exists 
only one such a point, s' , that is the solution to the equation 



s'h = Xz^M^s')). 



(4.28) 
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Specifically, we have for the pressure P = [3 V Ul8f . Using ( |4.24j ) and ( ft.lSQ we 
obtain an expression for the chemical potential in terms of (3 and P: 

where 

M 2 = m 2 + 2maP. (4.30) 



The standard thermodynamical relation |T6| 

dfi = -sdT + vdP, (4.31) 

where s = S/N and v = A/N , allows one to obtain the equation of state 

, d ± 1 | rna Kojm 

dP (3P M Ki((3M) ' y } 

It is convenient to represent this equation in the form: 

v=-£p + (r5(p,P), (4.33) 

where 

ruMPM} 

From the inequality 

K n (x) > K m (x), if n > m, (4.35) 

it follows that 

0<6(P,P)<1. (4.36) 

In the nonrelativistic approximation, when /3m — > oo, M — > m, the asymptotical 
expression ( p.llQ shows that 8 — > 1, and we obtain a well-known Tonks equation 
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(^) 

In the ultrarelativistic limit, when /3m — > 0, we have 5^0 and arrive at the 
state equation of the ideal gas. In the case of small a (that corresponds to the linear 
approximation in the interaction) we can replace function 5 ((3, P) by its value 5q((3) 
at a — > 0: 

% = (4-38) 

This quantity gets an interesting interpretation after accounting that the mean value 
of p n over the free-particle distribution defined by the Hamiltonian Q4.2p is given by 



(p k ) = / dpp k e-> JHoW / / dpe- pii °W =m k n+ Z ' ' . (4.39) 

Ki((3m) 



o 
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Then 5 (P) 

we have for free particles p 



m(p 1 )o- Next, we can observe that in the front form of dynamics 

-VT 



777.7, where 7 



2v corresponds to the Lorentz 



radical yl — v 2 . Therefore, in the linear approximation the state equation Q4.33|) 
can be considered as a result of taking into account the well-known Lorentz spatial 
contraction o — > cr(7 _1 )o in the nonrelativistic Tonks equation ( |4.37| ). 
From the equation of state ( |4.32|) we get immediately 



dv 
dP 



PP 2 



1 - 



Kim) 

K 2 ((3M)\ ' 



(4.40) 



Then inequality (|4.35|) shows that the condition of thermodynamical stability 



dv \ 

3pL <0 



(4.41) 



is valid for equation (|4.32| ) identically. Thus, the system does not have any phase 
transition. 

In a similar manner we may obtain explicit expressions for other thermodynam- 
ical functions. For example, for entropy we have 



df 



k0 



,dfi 

8(3 



k < In 



2MK 1 ((3M) 



h(3P 



the energy density is determined by 

[i + Ts — Pv 



u 



1 m 2 + maP K ([3M) 
Jj + M KxifiM) 



(4.42) 



(4.43) 



Specific heat cp can be directly obtained from ( f4.42| ) giving 



cp 



2 + (f3M) 2 



dp J 

K (PM) 2 

K^pM) 2 pMKtiPMy 



K (PM) 



(4.44) 



Therefore, we have an exactly solvable example of a nontrivial particle interaction 
in the relativistic statistical mechanics. 



5 Weakly-relativistic corrections to the thermo- 
dynamics of an interacting particle system 

Here we shall consider a system of N point-like particles with pairwise interactions at 
the first post-Newtonian approximation. Let the nonrelativistic interaction potential 
has the form: 

^ (0) = Uab; u ab = u(r a b), = \r ab \ = |x a -x b |. (5.1) 

a<b 
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In that case the general form of the first post-Newtonian Lagrangian function is 
given by 0: 

L = y- ( m a v l , rn a v£ \ _ (0) 
^ V 2 8c 2 / 

+ 77^- ( v a • \ b u ab - (r ab ■ v a )(r ab • v 6 ) — J + cT 2 $ (5.2) 

2 C ^ ^ j-^ y rab dr ab J 

a<o ' 

with an arbitrary Galilei-invariant function We choose this function in the form 
which is determined in terms of the nonrelativistic potential u(r): 



$ = - ^ ( Avl b u ab + B(r ab ■ v ab ) 2 — J , v afe = v a - v 6 , 

a<o 



(5.3) 



with two arbitrary numerical coefficients, A and B. This structure of the post- 
Newtonian Lagrangian is sufficient to cover a wide class of interactions including 
those following from various field-theoretical considerations. For example, the values 

A = a 2 - 1, B = (5.4) 

correspond to the interactions mediated by linear relativistic fields of spin a. Partic- 
ularly, cr = l and u(r) ~ r _1 give the famous Darwin's Lagrangian for electromag- 



netic interactions. It is remarkable that expression (|5.2|) under ( |5.3D was derived by 
Breit [^0| as far back as 1937 with the aid of simple symmetry treatment. 

For the system of identical particles the Hamiltonian which follows from the 
above Lagrangian has the form: 

H = # (0) + c- 2 # (1) +0(c- 4 ), (5.5) 

where is a nonrelativistic Hamiltonian, 

H (o) = yA + U (o) (5 . 6) 

a 

and 



a<6 



- [(1 + 2B)(r ab ■ p a )(r ab ■ p b ) - B[(r ab ■ p a ) 2 + (r ab ■ p b ) 2 ]} -1^1 . (5.7) 

Tab dr ab J 

The canonical variables (x a , p a ) are connected with the Lagrangian (x a , v a ) by the 
standard Legendre transformation 

Pa = (5.8) 

<9v„ 
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considered to the order c 2 . Inserting Hamiltonian (|5.5|) into expression ( |3.3|) for 
the classical partition function and expanding it to the c~ 2 terms, we get 



TV > 



where Z° N is a nonrelativistic partition function 



7 _ y(0)id n 



Q being a nonrelativistic configuration integral 

Q = V- N J e-^ C0) (1 - /3cr 2 if «) J] d 3 x a . 



(5.9) 



(5.10) 



(5.1i; 



The straightforward calculation gives first-order corrections to the nonrelativistic 
partition function (|5.10|) : 



7(1) _ 7(0) # . 
ZiV " ZiV /3m' 



7 _ 7(0) 



1 + 



1? 



/3mc 2 



(5.12) 



where i? is defined in terms of the configuration integral ( |5.11| ), 



(5.13) 



In the absence of interaction, when Q = 1, the obtained expression agrees with 
Eq. ([TT3D. 

All thermodynamical properties of the system may be deduced from the free 
energy: 

F = -p-HnZ = F^-— — 2 , (5.14) 
where the corresponding nonrelativistic expression is given by 



fW = -/3- 1 1b 



(5.15) 



Eq. ( |5.14j ) can be rewritten in the form, which gives a first-order correction in terms 
of the nonrelativistic free energy: 



p = f (0) - 



(3mc 2 



5 3A 



N 



-- — + B)- + A[ F<°> + (3 



dp ) 



BV 



«9K 



. (5.16) 



This formula may be also useful in the obtaining of weakly-relativistic corrections 
to various phenomenological nonrelativistic results. 

Let us consider corrections to the equation of state. In the nonrelativistic limit 
it has the form: 

P = -^ r =9 {0) (^V). (5.17) 
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Then (|5.16 ) gives a weakly-relativistic equation of state in the form 



j3mc 2 



(A + B)g n + A ^ + BV ^ 



(5.18) 



which is determined by the nonrelativistic form ( 5.1 7|) and the two constants, A and 
B. 

As an example we can consider a well known Van der Waals equation 



P 



N 



N 2 a 



p(y - Nb) v 2 



(5.19) 



with a pair of phenomenological constants, a and b. The corresponding weakly- 
relativistic equation is 



P 



N 



P(V - Nb) 



1 + 



3NbB 



f3mc 2 (y -Nb) 



N 2 a 
V 2 



1 + 



3(B - A) 

[3mc 2 



(5.20) 



Up to the first order in c 2 it can be presented in the nonrelativistic form ( |5.19|) 
with constants a and b being replaced by the linear functions on temperature: 



a \— >a' — a\l+3 



B- A 

(3mc 2 



b h+ b' = b 1 + 



3B 

f3mc 2 



(5.21) 



Corrections to other thermodynamical functions and a wider discussion can be found 



in Ref. 15 
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